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ABSTRACT 

The degradation of eukaryotic mRNAs can be initiated by deadenylation, decapping, or endonuclease cleavage. This is followed by 
5'-3' degradation by homologs of Xrn1, and/or 3'-5' degradation by the exosome. We previously reported that, in African 
trypanosome Trypanosoma brucei, most mRNAs are deadenylated prior to degradation, and that depletion of the major 5'-3' 
exoribonuclease XRNA preferentially stabilizes unstable mRNAs. We now show that depletion of either CAF1 or CNOT10, two 
components of the principal deadenylation complex, strongly inhibits degradation of most mRNAs. RNAi targeting another 
deadenylase, PAN2, or RRP45, a core component of the exosome, preferentially stabilized mRNAs with intermediate half-lives. 
RRP45 depletion resulted in a 5' bias of mRNA sequences, suggesting action by a distributive 3'-5' exoribonuclease. Results 
suggested that the exosome is involved in the processing of trypanosome snoRNAs. There was no correlation between effects 
on half-lives and on mRNA abundance. 
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INTRODUCTION 

The abundance of an mRNA in a cell is determined not only 
by its rate of synthesis and processing, but also by its decay 
rate. There are three ways in which degradation of an mRNA 
can begin: by removal of the 3' poly(A) tail or the 5' cap struc- 
ture (Li and Kiledjian 2010), or by endonucleolytic digestion 
(Schoenberg 2011). For most eukaryotic mRNAs that have 
been studied, constitutive degradation is initiated by deadeny- 
lation ( Chen and Shyu 2011). Three different types of cytosol- 
ic deadenylase are known. Poly(A) ribonuclease (PARN) is 
involved mainly in specialized processes (Chen and Shyu 
2011); PAN2 deadenylase, in a complex with PAN3, has a 
minor role, perhaps in initial shortening of poly(A) tails 
(Yamashita et al. 2005). The major activity resides in the 
Ccr4/Cafl/Not complex (Collart and Panasenko 2012), which 
consists of the central scaffold protein Notl, Not2, Not3 and/ 
or Not5, Caf40, either Cafl30 or CNOT10, the exoribonu- 
clease Cafl, and, in Opisthokonts, an additional exoribonu- 
clease, Ccr4 (Collart and Panasenko 2012). Shortening of 
the poly(A) tail to <20 nt leads to loss of poly(A) binding 
protein, which is followed by decapping. Some mRNAs are 
subjected to deadenylation-independent decapping, or to spe- 
cific endonuclease cleavage. In each case, the mRNA is then 
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exposed to 5' -3 degradation by the Xrnl exonuclease, or 3'- 
5' degradation by the exosome (Chen and Shyu 201 1). 

Models of eukaryotic mRNA decay were established using 
a very small selection of specific transcript sequences but 
have been followed by transcriptome-wide analyses. Deletions 
of Ccr4 and Pan2 in yeast caused stabilization of mRNAs en- 
coding ribosomal proteins and ribosome biogenesis factors 
(Grigull et al. 2004): These mRNAs are among the least stable 
in yeast. More recently, Saccharomyces cerevisiae mRNA decay 
rates were measured by 4-thio-uracil pulse-chase labeling, and 
a median half-life of 11-12 min was measured (Miller et al. 
2011). Interestingly, cells expressing a mutant RNA polymer- 
ase II with slow elongation compensated by decreasing mRNA 
degradation; conversely, cells lacking Ccr4p or Caflp had 
globally decreased mRNA degradation and compensated via 
decreased mRNA synthesis (Sun et al. 2012). In the latter 
case, deletion of CAF1 had the stronger effect, since it abolishes 
deadenylation by the Not complex: The average half-life in- 
creased over sixfold, while synthesis decreased 2.6-fold (Sun 
etal. 2012). 

The exosome core has nine subunits: a central hexamer 
surmounted by three SI -domain proteins. Exosome core 
components are essential, since the exosome is required for 
processing of many stable RNAs including rRNA, in addition 
to mRNA. Exosome RNase activity is contributed by associ- 
ated 3'-5' exoribonucleases: Rrp6p (mainly restricted to the 
nucleus in mammals and yeast) and Rrp44p (mainly in the 
cytoplasm). Shutting off expression of a core exosome com- 
ponent in yeast caused increases in mRNAs for ribosome 
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biogenesis, processing of 20S RNA, and nucleo-cytoplasmic 
transport (Mnaimneh et al. 2004). Reduction of a different 
core component, however, resulted in a decrease in mRNAs 
encoding ribosomal components, which was attributed to 
growth arrest (Houalla et al. 2006). Depletion of an exosome 
subunit in Drosophila cells caused increases in steady-state 
levels of ~1000 different mRNAs (Kiss and Andrulis 2010); 
mRNAs involved in the cell cycle and development were 
preferentially affected, as were mRNAs with unusually long 
untranslated regions. Only 60 mRNAs were decreased. In- 
ducible RNAi against Arabidopsis exosome subunits revealed 
increases in ~250 specific mRNAs, as well as various novel 
RNAs, and RNAs with unexpected 3' extensions (Chekanova 
et al. 2007). The exosome is involved in the processing of 
snRNAs and snoRNAs; the absence of Rrp6p or other exo- 
some subunits results in an accumulation of precursors, in- 
cluding polyadenylated degradation intermediates, in both 
yeast (Wyers et al. 2005; Gudipati et al. 2012) and Arabidopis 
(Chekanova et al. 2007). Increases in mRNAs located down- 
stream from snoRNA loci were also observed (Houalla et al. 
2006). 

The transcription of protein-coding genes in Kinetoplastid 
protists by RNA polymerase II is polycistronic (Palenchar and 
Bellofatto 2006; Daniels et al. 2010). Individual mRNAs are 
excised by 5' trans splicing of a "spliced leader" (SL) and 3' 
polyadenylation (Michaeli 201 1). Since transcription of indi- 
vidual genes cannot be controlled, mRNA turnover plays an 
unusually dominant role in determining mRNA abundance 
(Ouellette and Papadopoulou 2009; Fernandez-Moya and 
Estevez 2010; Kramer and Carrington 2011; Manful et al. 
201 1). For most mRNAs, removal of the poly(A) tail precedes 
degradation of the remainder of the transcript (Manful et al. 
2011). A typical mRNA degradation machinery is present: 
an exosome (Estevez et al. 2001, 2003); a scavenger decapping 
enzyme (Banerjee et al. 2009); an Xrnl homolog, XRNA 
(Li et al. 2006); PAN2/PAN3 (Schwede et al. 2009); PARN 
(Milone et al. 2004; Utter et al. 201 1); a CAF1/NOT complex 
composed of CAF1, CAF40, NOT1, NOT2, NOT3/5 (Sch- 
wede et al. 2008); a CNOT11 homolog (Schwede et al. 2008; 
Mauxion et al. 2013); and CNOT10 (Farber et al. 2013). 
Only a Dcp2-like decapping enzyme is missing (Schwede 
et al. 2009). RRP45 is an essential structural component of 
the exosome; its depletion results in complex disassembly 
and instability of other exosome subunits (Estevez et al. 
2001, 2003). Similarly, CNOT10 is required for the attach- 
ment of the deadenylase CAF1 to the CAF1/NOT complex 
(Farber et al. 2013). 

The median mRNA half- life in bloodstream- form trypano- 
somes (the form that grows in mammals) is ~ 15-20 min 
(Manful et al. 2011). Depletion of XRNA (Li et al. 2006), 
PAN2 (Schwede et al. 2009), or the exosome (Haile et al. 
2003) delayed the degradation of the highly unstable develop- 
mentally regulated mRNA encoding EP procyclin, or reporter 
mRNAs bearing its 3' UTR; ACT (actinA) mRNA, which has a 
half-life at the median, was also stabilized by RNAi targeting 



PAN2, and the developmentally regulated unstable mRNA 
encoding pyruvate phosphate decarboxylase was stabilized 
by RRP45 RNAi. In contrast, depletion of CAF1 (Schwede 
et al. 2008, 2009) or CNOT10 (Farber et al. 2013) affected 
all mRNAs examined. Detailed studies of EP mRNA degrada- 
tion kinetics provided clear evidence for a very rapid, deadeny- 
lation-independent, and XRNA-dependent 5'-3' pathway, 
combined with a slower deadenylation-dependent compo- 
nent (Irmer and Clayton 2001; Haile et al. 2008; Schwede 
et al. 2009). In the related Kinetoplastid parasite Leishmania 
there is also evidence for this type of biphasic degradation 
(Haile et al. 2008), and some developmentally unstable 
mRNAs are subject to endonuclease cleavage in the 3' UTR 
(Muller et al. 2010). So far, endonucleolytic pathways have 
not been found in trypanosomes. 

Transcriptome-wide analysis of trypanosome mRNA half- 
lives in cells after XRNA RNAi demonstrated that XRNA 
depletion inhibited the degradation of unstable mRNAs, 
while some of the more stable mRNAs had shorter half- 
lives (Manful et al. 2011). We have now extended this to all 
the known factors involved in degradation from the 3' end: 
PAN2, the exosome, and components of the CAF1/NOT com- 
plex. Our results confirm that the CAF1/NOT complex is 
central to RNA turnover in trypanosomes, while depletion 
of the exosome and PAN2 complexes mainly delayed degrada- 
tion of mRNAs with median half-lives. 

RESULTS 

The rate of bulk mRNA degradation in wild-type cells 

The aim of this work was to assess the transcriptome-wide role 
of 3' pathways in mRNA degradation, using RNA-seq. We 
first, however, show results for bulk mRNA. Given the unusual 
mode of trypanosome mRNA synthesis and processing, to in- 
hibit creation of new mRNAs, we need to inhibit not only tran- 
scription, but also processing of preexisting precursors. To 
inhibit processing, we therefore first added Sinefungin, which 
inhibits methylation of the 5' cap of the SLRNA used as a donor 
for the 5' 39-nt spliced leader (SL). Unmethylated SLRNA is 
not a substrate for trans splicing (Ullu and Tschudi 1991). 
Five minutes later, when we expected reduced cap methyla- 
tion, we added Actinomycin D (Colasante et al. 2007). 

Since all trypanosome mRNAs carry the SL at the 5' end, the 
abundance of the SL can serve as a proxy for the abundance of 
total mRNA, independently of the presence of a poly(A) tail. 
To determine the average half-life of the bulk mRNA popula- 
tion, we therefore hybridized Northern blots with an SL probe. 
As previously shown (De Lange et al. 1984), a smear was ob- 
tained, with more prominent areas at the positions of abun- 
dant mRNAs. Figure 1 shows results for samples that were 
later subjected to RNA-seq, while Supplemental Figure SI 
shows more detailed time courses. The SLRNA (arrow, Fig. 
1 ) is the 5' substrate for trans splicing, which consists of the 
5' SL and the SL intron. The SLRNA disappeared rapidly after 
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FIGURE 1 . Spliced leader hybridization to Northern blots of the RNA 
samples sent for sequencing. RNA from different cells was prepared and 
subjected to Northern blotting with a spliced leader RNA 39mer oligo- 
nucleotide probe. Three similar samples were tested for each line; one 
typical lane is shown. A second set of lanes from wild-type (WT) and 
CNOT10 RNAi was shown in Farber et al. (2013). The arrow points 
to the spliced leader precursor RNA, and the asterisk marks signal 
from the 7SL (signal recognition particle) RNA left from a previous hy- 
bridization. The mRNA signal was quantitated in the region from ~400 
nt to 6 kb. The panel below shows the methylene blue signal from rRNA 
on the corresponding blot segments. The hybridization signals for indi- 
vidual RNAi lines differ because RNA from each cell line was run on a 
different gel. For further comparisons see Supplemental Figure SI. 



Actinomycin D addition to wild- type cells; SLRNA with meth- 
ylated caps could be used in splicing, while the remainder was 
presumably degraded (Fig. 1; Supplemental Fig. SI). To assess 
mRNA levels, the density of the smear >400 nt was measured, 
then relative amounts were calculated using rRNA or the 7SL 
(signal recognition particle RNA) as a standard. The results for 
RNA preparations that were subsequently used for RNA-seq 
(Fig. 1) are listed in Table 1. The proportion of SL signal 
that remained after 30-min transcription inhibition in wild- 
type cells was 50%-60% (Table 1), which was approximately 
in agreement with our previous estimates based on the abun- 



dances of individual mRNAs (Haanstra et al. 2008; Manful 
etal. 2011). 

Supplemental Figure SI A shows the time course of disap- 
pearance of SL signal over a 2-h time course. Three inde- 
pendent experiments yielded half-lives between 44 and 46 
min, with 54%-56% of SL signal remaining after 30 min 
(Supplemental Fig. S2). 

The rate of bulk mRNA degradation in cell lines 
deficient in 3 -5' RNA degradation pathways 

To analyze the roles of deadenylation and the exosome in 
global mRNA decay, we used cell lines in which the mRNAs 
encoding CAF1 (Tb927.6.600), CNOT10 (Tb927.10.8720), 
PAN2 (Tb927.6.1670), or RRP45 (Tb927.6.670, an essen- 
tial core component of the exosome) were depleted. For the 
RNA-seq analyses, we used existing cell lines (Estevez et al. 
2001; Schwede et al. 2008, 2009; Farber et al. 2013). All lines 
were checked for the effectiveness of RNAi, by Northern blot- 
ting and by ensuring that growth was inhibited within 48 h 
of RNAi induction (data not shown) , and later, more accurate- 
ly, by analysis of the RNA-seq results (Table 1 ) . For all exper- 
iments, RNAi was induced for 24 h; this time point precedes 
the onset of growth inhibition. Analysis of the amounts of 
SL-bearing mRNA that remained after 30-min Actinomycin 
D indicated that the bulk mRNA decay rates in cells with 
PAN2 RNAi differed little from wild type, whereas for CAF1, 
mRNA degradation was almost completely inhibited (Fig. 
1; Table 1). The strong inhibition of mRNA decay after 
CNOT10 RNAi was previously reported (Farber et al. 2013) 
and is illustrated again in Table 1, Figure 1, and Supplemental 
Figure SI. Bulk mRNA degradation in RRP45 RNAi lines is 
discussed in the next section. 

Results for 2-h time courses of SL-mRNA degradation are 
shown in Supplemental Figure SI and quantitation is in 
Supplemental Figure S2. These experiments were done later 
than those in Figure 1 . In part, new cell lines had to be made 



TABLE 1. Characteristics of mRNA degradation in exoribonuclease mutants 



Cell line 


WT 


RRP45 


PAN2 


CAF1 


CNOT10 


XRNA 


Average fold target mRNA decrease 3 




3x 


11x 


2x 


82 x 




SL mRNA at 30 min b 


55 ± 22% 


ND 


60 ± 2% 


93 ± 1 9% 


1 02 ± 1 4% 




TUB mRNA at 30 min b 




76 ± 5% 


ND 


-100% 


~1 00% 




HISH4 mRNA at 30 min b 




66 ± 6% 


ND 


-100% 


~1 00% 




No. of mRNAs with half-life >1 20 min a 


146 


199 


506 


1616 


4471 


16 


Median bulk mRNA half-life (min) a 


12 


14 


15 


32 


>240 


18 


Arithmetic mean bulk mRNA half-life (min) a 


21 


25 


39 


90 


>200 


21 



Results for wild-type and XRNA-depleted cells are from Manful et al. (201 1). 

a Results from RNA-seq; for the calculations only the 5326 "unique genes" of Siegel et al. (2010), and only ORFs with an RPM of at least 5 at 
steady state were included. 

b These are the results for 0- and 30-min time points, measured by Northern blotting, for the samples that were prepared for sequencing. The 
SLRNA signals were normalized to the 7SL RNA signal but were also checked using the methylene blue rRNA stain. If a standard deviation is 
shown, the data are for three samples prepared for RNA-seq, of which two were actually sequenced; If no standard deviation is shown (CAF1 
and CNOT10, TUB and HISH4), only two samples are included. No results are shown for XRNA since the data used for comparisons are 
from an earlier paper (Manful et al. 201 1). 
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because the effectiveness of the RNAi in 
the stored cell lines had diminished. This 
means that the results for RNA-seq (Fig. 
1; Table 1) and for the 2-h time courses 
can be compared only qualitatively. In 
the 2-h experimental series, RNAi tar- 
geting XRNA, PAN2, and RRP45 had 
no consistent effect on bulk mRNA 
turnover (Supplemental Figs. S1B-D, 
S2), whereas, once again, depletion of ei- 
ther CAF1 or CNOT10 strongly delayed 
mRNA degradation (Supplemental Figs. 
S1E.F, S2). RNAi targeting CAF1 was, in 
these experiments, somewhat more ef- 
fective than RNAi targeting CNOT10 
in inhibition of mRNA decay; this can 
be attributed to the use of new cell lines 
with different RNAi penetrance. 

Depletion of RRP45 results in a 5' 
read bias 
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We used RNA-seq to examine transcrip- 
tomes at steady state and 30 min after 
inhibition of mRNA maturation and 
synthesis with Sinefungin and Actino- 
mycin D, respectively. From the three 
replicates used to estimate the disappear- 
ance of the SL signal (Table 1), two were 
sent for sequencing. All raw data (read 
counts) are in Supplemental Table 1 
and the numbers of aligned reads are 
shown in Supplemental Figure S3. We 
intended to use the SL hybridization sig- 
nal (examples in Fig. 1 ) to normalize our 
RNA-seq data, but this is only valid if 
the 5' end is representative of the whole 
mRNA. We therefore examined the dis- 
tributions of read counts across the mRNAs both before and 
after addition of inhibitors. For this purpose, we extracted 
the most commonly used 5' splice and 3' polyadenylation sites 
from our data. Next, for every mRNA that had mapped ex- 
tremities, we normalized the read density results to a stan- 
dardized length of 1 kb. This enabled us to plot the average 
read densities for all mRNAs irrespective of length. Results 
are shown in Figure 2. For untreated cells there was a mild 
decrease in density toward the 3' end (Fig. 2). Sometimes, 
this decrease already started within the mRNA coding region 
(data not shown). In the wild-type cells, after 30-min tran- 
scription inhibition, the read density pattern had hardly 
changed; the mRNA population at 30 min therefore presum- 
ably consisted mainly of intact mRNAs (Fig. 2). 

Next, we compared the distributions for the RNAi cell lines, 
incorporating also our previous results from XRNA RNAi 
(Manful et al. 201 1). The overall distributions, both at steady 
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FIGURE 2. Depletion of RRP45 results in a 5' read bias. Predominant polyadenylation and splice 
acceptor sites for as many mRNAs as possible were extracted from all our data sets. For each 
unique ORF, the read density was plotted across the entire predicted mRNA, then normalized 
to an mRNA length of 1000 nt. The average read density across all mRNAs was then plotted. 
Results for RNAi lines are compared with wild type, as indicated on the individual panels (A-E). 



state and after 30-min Actinomycin D, were in general similar 
to wild type (Fig. 2A-D). However, for cells with RRP45 RNAi 
after 30-min Actinomycin D, there was a clear readbias toward 
the 5' end (Fig. 2E). This suggests that the cells contained 
mRNAs that were in various stages of 3'-5' degradation. For 
these cells, therefore, we could not use SL signals in order to 
normalize the data, and relied on a tubulin instead (see below) . 

Inhibition of 3 -5' pathways preferentially affects 
more stable mRNAs 

We next used the transcriptome data to obtain estimates of 
mRNA half-lives for all ORFs, assuming simple exponential 
decay (Supplemental Table S2) . Since the SL results were prob- 
ably biased in the RRP45 samples, we normalized those based 
on degradation of the a tubulin mRNA (Table 1), which was 
very similar to wild type (data not shown). The remainder 
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were normalized using the SL Northern blot measurements 
from the original RNAs and additional replicates that were 
made at the same time (Fig. 1 ; Table 1 ) . We selected all unique 
ORFs, and added up the total read counts obtained from them. 
From these, we calculated reads per million (RPM) for all 
ORFs. For example, for PAN2, in which 60% of total mRNA 
remains after 30 min (Table 1), all RPMs from the 30-min 
time points were multiplied by 0.6. The 30-min RPM for 
CAF1 and CNOT10 were not changed since no significant 
degradation of total mRNA had occurred after 30 min. 
Results for the duplicates were reproducible (log scale correla- 
tion coefficients exceeding 0.85), with most variation being 
concentrated at low read densities (Supplemental Fig. S3). 
For most subsequent calculations, we therefore considered 
only mRNAs represented by at least 5 RPM at time 0 in the 
wild-type sample. The use of a single time point gives only a 
very crude estimate of the half-life: Complex kinetics 
(Deneke et al. 20 1 3 ) are not detected, and since we are extrap- 
olating, the more the half-life differs from 30 min, the greater 
the error is likely to be. For the figures and tables the maximum 
half-life was arbitrarily set to 240 min. 

Despite the limitations of our approach, some conclusions 
can be drawn. The clearest result was that depletion of PAN2 
(Fig. 3A) or RRP45 (Fig. 3B) had much less effect on the 
half-life distribution than depletion of CAF1 (Fig. 3C) or 
CNOT10 (Fig. 3D). Impairing CAF1/NOT complex function 
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FIGURE 3. Effect of RNAi on the half-lives of trypanosome mRNAs. 
The half-life of the mRNA from each unique ORF in cells with RNAi 
is plotted against the half-life that was previously measured for wild- 
type cells, on a log scale. The dashed diagonal indicates perfect correla- 
tion. All half-lives >240 min were set to 240 min, and all ORFs with a 
steady-state RPM of <5 for any of the cell lines were excluded from 
the data set. RNAi was as follows: (A) PAN2; (B) RRP45; (C) CAF1; 
(D) CNOT10. 
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FIGURE 4. Effect of RRP45, PAN2, and CAF1/NOT complex depletion 
on mRNA decay. (A) Effects on the overall half-life distribution. Each 
unique ORF was placed in a half-life category as indicated in the key 
on the right. The numbers of ORFs in each category are displayed. (B) 
The half-life for each unique ORF after RNAi was divided by the half- 
life in wild type. The ORFs were then organized into categories depend- 
ing on the wild-type half-life, and the average change in half-life after 
RNAi was calculated. The differences between the means of the catego- 
ries for all RNAi lines were significant (one-way ANO VA on linearized 
dataP-value <1 x 10" 16 ). 



delayed degradation of almost all mRNAs, with most now 
having half-lives >60 min. Thus it was clear that, for the vast 
majority of mRNAs, deadenylation by the CAF1/NOT com- 
plex is a critical event in the onset of mRNA decay. CNOT10 
is required for attachment of CAF1 to the remaining NOT 
subunits and for CAF1 deadenylation activity in vivo (Farber 
et al. 2013). CNOT10 RNAi had a stronger effect on mRNA 
decay than CAF1 RNAi, but CNOTW RNA was also more 
effectively depleted (>10x reduction as opposed to only 2x 
for CAF1). Moreover, CAF1 protein is present in excess 
(Farber et al. 2013). In subsequent discussions, we will there- 
fore treat the two sets of results together, assuming that the 
difference between the two is only due to variations in the ef- 
fectiveness of the RNAi. 

We next looked at the distributions of half-lives. To do this, 
each unique ORF was placed into a half-life category, then the 
number of ORFs in each category was counted (Fig. 4A). In 
wild-type cells, we had previously found that most mRNAs 
had half-lives <20 min, with very few >1 h (Fig. 4A; Manful 
et al. 201 1). Depletion of PAN2 and RRP45 had little effect on 
the overall distribution, whereas disabling the CAF1/NOT 
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complex shifted many mRNAs into the longer-lived categories 
(Fig. 4A). To see whether mRNAs with particular characteris- 
tics were specifically affected by inhibition of different 3'-5' 
degradation pathways, we first calculated for each ORF the 
ratio between the wild-type half-life and the new half-life. 
We then grouped mRNAs according to their half-lives in 
wild-type cells, and plotted the average ratio for each group. 
Note that, because of our arbitrary maximum half- life of 240 
min, an mRNA with a wild- type half-life of 120 min cannot 
increase more than twofold in the tables, so effects on more 
stable mRNAs may be underestimated. We had previously 
found that, after XRNA depletion, mRNAs with short half- 
lives were specifically stabilized, whereas mRNAs with longer 
half-lives were destabilized (Manful et al. 2011); this is again 
illustrated in Figure 4B. The effects of PAN2 and RRP45 deple- 
tion were different: mRNAs with intermediate half-lives were 
stabilized, but again, for RRP45, the long-half-life mRNAs 
appeared to be less stable (Fig. 4B). CAF1 and CNOT10 de- 
pletion also affected the least stable mRNAs less than those 
of intermediate stability. The significance of these different ef- 
fects for all cell lines was confirmed by one-way ANO VA. Since 
many mRNAs with wild- type half-lives >20 min had no mea- 
surable degradation in these two mutants, no data for these 
are plotted because the fold effect could not be calculated. 

Changes in mRNA decay are not reflected 
in steady-state abundance 

We next examined the effects of RNAi on the steady- state 
abundances of mRNAs. There was no effect on bulk mRNA 
abundance, as judged by the SL signal (Supplemental Fig. 
1G). We had previously seen that, after XRNA depletion, un- 
stable mRNAs became slightly more abundant than in wild- 
type (Manful et al. 2011), although the correlation was very 
weak. The drastic inhibition of mRNA decay seen after inhibi- 
tion of the CAF1/NOT complex would be expected to have 
far greater effects. When, however, we examined the abun- 
dances of mRNAs after RNAi against CNOT10, RRP45, or 
PAN2 there was no correlation at all between the changes in 
half-life and those in mRNA abundance (Fig. 5; Table 2; 
Supplemental Table S3). For RRP45 and PAN2 knockdowns, 
there were many ORFs showing altered mRNA abundance de- 
spite little change in half-life. To examine this further, for the 
CAF1/NOT complex, we looked at mRNAs that had half-lives 
of <240 min in both wild- type and CAF1 RNAi cells. The pres- 
ence of a significant expression change for an ORF (P-adjust- 
ed- value <0.01) was not predictive of an alteration in mRNA 
decay rate. 

To check whether the discrepancy between change in abun- 
dance and change in half-life was in any way dependent on 
the position of genes in a transcription unit, we divided the 
change in abundance by the change in half-life for each unique 
gene. We then, for chromosomes 1,2, and 3, plotted the results 
in the correct gene order. The distribution of effects appeared 
to be random, both between transcription units and with 
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Log, 0 half-life ralio (CNOT10 RNAiWT) 




Log 10 half-life ratio (PAN2 RNAiWT) 




-15 

log,o half-life ratio (RRP45 RNAiAVT) 



FIGURE 5. Changes in mRNA half-lives after RNAi do not lead to cor- 
responding changes in mRNA abundance. The mRNA half-life after 
RNAi (in minutes) was divided by that in wild-type cells and plotted 
on the x-axis (log 10 ); the mRNA abundance after RNAi (average RPM 
reads) was divided by that in wild-type cells and plotted on the y-axis 
(log 10 ). mRNAs with half-lives >240 min were excluded. Trend lines 
were plotted in Excel. (A) CNOT10 RNAi; (B) PAN2 RNAi; (C) 
RRP45 RNAi. 

respect to position of a gene relative to the transcription start 
site. 

RNAi targeting 3 -5' decay preferentially affects mRNAs 
encoding proteins with specific functional classes 

We next looked to see whether the steady-state levels, or half- 
lives, of mRNAs encoding particular functional classes of pro- 
teins were differentially affected by inhibition of 3'-5' mRNA 
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TABLE 2. Statistics for genes showing changed expression at steady state 





Number of 


Mean fold 


change in 


Mean WT half-life 




differentially expressed 


half-life, for 


genes with 


(min) for ; 


^enes with 




8 


enes 


mRNA that was: 


mRNA that was: 


RNAi target 


increased 


decreased 


increased 


decreased 


increased 


decreased 


CNOT10 


96 a 


305 a 


6.0 


3.2 


45 


55 


CAF1 


714 


708 


3.0 


9.0 


32 


30 


PAN2 


439 


761 


1.5 


1.3 


34 


24 


RRP45 


1049 


1224 


1.0 


2.8 


25 


30 



We included only unique genes (according to Siegel et al. [2010] with a few manual 
additions), with a minimum of 5 RPM at steady state and 1 RPM after 30-min Actinomycin 
D. The mRNA abundance was judged to be increased or decreased when it was statistically 
significant. 

a The low number of differentially expressed genes for CNOT1 0 is a technical artifact, as for 
many genes the difference between replicates was too big to allow a clear statement to be 
made. 



snoRNAs themselves are much too small 
to be reliably retained in our data set, but 
precursors of >250 nt could have been 
detected. It was therefore interesting 
that there were far more snoRNA reads 
from the RRP45 RNAi cells than from 
either wild type or any of the other 
mutants (Supplemental Table S3; Sup- 
plemental Figs. S4, S5). Examination 
of the raw data (Supplemental Figs. 
S4, S5) confirmed that RNA sequence 
reads from snoRNA clusters were spread 
across both the mature snoRNA and 
the neighboring regions. This strongly 
suggests that the exosome is involved 
in snoRNA processing in trypanosomes, 
as was previously demonstrated in other 
species (see Introduction). 



decay. To classify trypanosome genes, we supplemented the 
current annotations from GeneDB and TritrypDB with 
KEGG designations and a limited number of publications de- 
scribing pathways and organellar or complex proteomes 
(Supplemental Table S2; Bringaud et al. 2006; Colasante et 
al. 2006; Michels et al. 2006; Panigrahi et al. 2008; Vertom- 
mena et al. 2008; Acestor et al. 2009; Droll et al. 2013). We 
had already seen that, in wild- type cells, mRNAs encoding ri- 
bosomal proteins and enzymes of glucose metabolism are sig- 
nificantly more stable than the average (Manful et al. 201 1). 

First we examined steady-state mRNA levels. In all of the 
mutants, ribosome-related genes (mostiy encoding ribosomal 
proteins) were significantly enriched in the group of mRNAs 
that were significantly increased in abundance (Fisher test P- 
value <1 x 10~ 8 ) (Supplemental Table S3). In contrast, cyto- 
skeleton-related genes were significantly enriched in the group 
showing significantly decreased mRNA abundance (Fisher 
testP-value <1 x 10" 8 ) (Supplemental Table S3). 

We now looked for functional classes that were enriched in 
the group of mRNAs that had an increase of at least twofold 
in their half-life (Fig. 6). For both RRP45 and PAN2 RNAi, 
mRNAs encoding cytoskeletal proteins were very signifi- 
cantly enriched in the stabilized category, and those encoding 
chaperones were stabilized by PAN2 RNAi. The most notably 
destabilized categories were mRNAs encoding enzymes of 
glucose and glycerol metabolism, and RNA-binding proteins 
(Fig. 6). Developmentally regulated mRNAs were not specif- 
ically stabilized by depletion of any of the 3' -5' decay pathway 
components (data not shown). No general enrichment anal- 
ysis was done for the NOT complex mutants since the major- 
ity of mRNAs were stabilized. 

RNAi targeting RRP45 results in increased 
read counts from snoRNA precursors 

Many trypanosome snoRNAs are processed from longer 
precursors (Doniger et al. 2010; Michaeli et al. 2012). The 



DISCUSSION 



We have shown here that a reduction in CAF1/NOT complex 
activity, effected through depletion of mRNA encoding either 
the CNOT10 or CAF1 subunit, inhibits the degradation of 
most trypanosome mRNAs. Depletion of either RRP45 or 
PAN2 affected much smaller subsets of mRNAs, which most- 
ly did not overlap with each other (Fig. 7). Changes in mRNA 
half-life did not generally correlate with alterations in mRNA 
abundance, showing that additional regulatory mechanisms 
must be present. 

In yeast with mutations in mRNA degradation, changes in 
mRNA decay were not reflected in mRNA steady-state abun- 
dance (Sun et al. 2012), and the suggested explanation was 
that feedback controls on transcription compensate for the 
lack of mRNA decay. The lack of correlation in trypanosomes 
is much more difficult to understand, because trypano- 
somes cannot alter the transcription rates of individual genes. 
Transcription by RNA polymerase II was clearly still active 
at the time we did our experiments; otherwise, a 30-min incu- 
bation with Actinomycin D would have had no effect at all on 
the relative abundances of different mRNAs and no mRNAs 
would have appeared to retain short half-lives. The most ob- 
vious possibility is that changes in trans splicing compensate 
for defects in mRNA degradation. 

RNAi targeting PAN2 and RRP45 accelerated the degrada- 
tion of quite a large number of mRNAs. These were enriched 
in categories that are required for growth, such as ATP gen- 
eration from glucose and glycerol, vesicular transport, nucle- 
otide metabolism, and translation. One possibility is that the 
effects were secondary to the onset of growth inhibition. It 
was notable that RNA-binding proteins were also particularly 
enriched in the destabilized category; perhaps some changes 
in mRNA stability or abundance reflected loss of regulatory 
RNA-binding proteins. 

The read densities at steady state were not uniform across 
the whole length of mRNAs. Instead, densities were roughly 
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FIGURE 6. The effects of RNAi against PAN2 and RRP45 on genes en- 
coding proteins of different functional classes. All unique ORFs were as- 
signed into categories according to the function of the encoded protein. 
Those with mRNA stabilization or destabilization of at least twofold 
were selected, and the relative enrichment of each functional class was 
then calculated. Rows represent gene classes; columns represent the sta- 
bilized or destabilized mRNA. Each rectangle is colored to represent the 
log P-values for the enrichment of the corresponding functional cate- 
gory, according to Fisher's test. The darker the pink color, the more sig- 
nificant the enrichment (see key). 



most frequently used poly(A) site will often be downstream 
from several other, less frequently used sites, so that the 
steady-state mRNA population is certain to contain some 
mRNAs with shorter 3' UTRs. 

After RNAi targeting PAN2, CAF1, or CNOT10, the distri- 
butions of sequence reads across RNA lengths were similar 
to those for wild type. This is the expected result if the RNAi 
causes accumulation of intact (polyadenylated) mRNAs. 
Depletion of XRNA did not affect the profile either for all 
genes (Fig. 2) or for those mRNAs whose degradation was 
clearly delayed by XRNA depletion (data not shown). This 
strongly suggests that trypanosome XRNA is a processive en- 
zyme. In contrast, RNAi targeting RRP45, which results in im- 
paired exosome activity, caused a change in the profile that 
suggested an accumulation of degradation intermediates con- 
taining the 5' ends of mRNAs. This could only happen if 
mRNAs had been deadenylated, then subjected to the action 
of limiting amounts of a distributive 3'-5' exoribonuclease. 
We do not know why these mRNAs were not then degraded 
by XRNA instead; interactions between the exosome and the 
decapping complex or XRNA cannot be ruled out. The distrib- 
utive 3'— 5' activity could come either from limiting amounts 
of the exosome or from a different 3'-5' exoribonuclease. 

Cytoplasmic mRNA that is targeted by Opisthokont exo- 
some is thought to be threaded through a hole in the exosome 
core, then degraded by the exoribonuclease Rrp44/Dis3, 
which is situated at the pore exit (Wang et al. 2007; Bonneau 
et al. 2009). Rrp44 is a processive enzyme (Lorentzen and 
Conti 2005; Liu et al. 2006; Dziembowski et al. 2007). So far, 
however, there is no evidence for an association between the 
trypanosome exosome and the trypanosome Rrp44 homolog 
(Estevez et al. 2001, 2003). Instead, trypanosome exosomes 
contain a homolog of Rrp6 (Estevez et al. 2001, 2003; Haile 
et al. 2007). In yeast and mammals, Rrp6 is associated pre- 
dominantly with nuclear exosomes and, interestingly, it is a 
distributive enzyme (e.g., Liu et al. 2007; Januszyk et al. 
2011). Evidence from cryo-electron microscopy and X-ray 
crystallography also so far suggests that RRP6 is located on 
one side of the exosome, rather than at the pore exit 
(Cristodero et al. 2008; Makino et al. 2013). It is therefore pos- 
sible that, in trypanosomes, exosomal degradation of mRNAs 
is indeed distributive. 



uniform across the first 2/3, then gradually declined toward 
the 3' end. This could be a technical artifact. Both gene length 
and nucleotide composition can affect the patterns of RNA 
fragmentation and fragment amplification, resulting in clear 
biases in RNA-seq results (see for example Dohm et al. 
2008; Sendler et al. 2011; Zheng et al. 2011). We have, how- 
ever, found no base bias in our results (Manful et al. 2011). 
Another contributing factor could be the fact that possible 
polyadenylation sites in trypanosome mRNAs can be distrib- 
uted over large distances — in extreme cases, several kilobases 
(Kolev et al. 2010; Nilsson et al. 2010; Siegel et al. 2010). The 




FIGURE 7. Some mRNAs are stabilized in all RNAi lines. The ORFs 
with at least twofold increased half-lives were selected for each RNAi 
line. The Venn diagram shows the overlaps between the different lines. 
The complete gene set is represented by the largest circle. 
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Our previous studies of 3'-5' mRNA decay in bloodstream- 
form trypanosomes focused on very few transcripts (Haile 
et al. 2003; Schwede et al. 2008, 2009; Farber et al. 2013). 
The data described here show that results from individual 
mRNAs cannot be generalized in order to make conclusions 
about all mRNAs with a particular half-life. Our results are, 
however, consistent with previous evidence that the 5' -3' 
and 3'-5' pathways can operate independently on a relatively 
unstable mRNA (Irmer and Clayton 2001; Haile et al. 2003; 
Li et al. 2006; Schwede et al. 2009). We already saw that 
many of these are stabilized after XRNA depletion (Manful 
et al. 2011). Our current finding that CAF1 and CNOT10 
depletion had less effect on degradation of the least unstable 
mRNAs than on those with intermediate stability (Fig. 4B) 
also fits the model. Most of the mRNAs that showed at least 
twofold increases in half-life after XRNA depletion (Manful 
et al. 2011) were affected by CAF1 or CNOT10 depletion 
but there was very little overlap between these and those that 
were stabilized after loss of RRP45 or PAN2 (Fig. 7). All 
mRNAs that were stabilized after PAN2 depletion, and the ma- 
jority of those that were stabilized after RRP45 depletion, were 
also dependent on the CAF 1 /NOT complex and there was a lit- 
tle under 50% overlap between these two gene sets (Fig. 7). 
This makes sense for the exosome, which can only act after 
poly (A) -binding protein has been removed from the 3' end; 
and it is consistent with other data suggesting that PAN2 alone 
does not normally effect complete deadenylation (Chen and 
Shyu2011). 

MATERIALS AND METHODS 
Cell lines and culture 

Trypanosoma brucei bloodstream forms expressing the tet repressor 
were cultivated as described previously (Alibu et al. 2004). Cell lines 
used had tetracycline-inducible RNAi targeting CAF1 (Schwede et al. 
2008), CNOT10 (Farber et al. 2013), PAN2 (Schwede et al. 2009), or 
RRP45 (Estevez et al. 2001). Cells were grown without selecting anti- 
biotics for 24 h, then RNAi was induced by addition of 1 00 ng/mL tet- 
racycline for 24 h. Tetracycline treatment alone does not affect the 
transcriptome (Manful et al. 20 1 1 ) . For analysis of RNA degradation, 
cells were treated with Sinefungin ( final concentration 2 pg/mL) for 5 
min, then Actinomycin D was added (final concentration 10 pg/mL) 
and cells were incubated for a further 30 min for the RNA-seq anal- 
ysis, or for additional times for the time courses. The time taken to 
centrifuge the cells (at room temperature) was included as part of 
the incubation time. 

RNA analysis 

Total RNA was extracted using peqGold Trifast (peqLab). For quality 
control RNA was run on a formaldehyde gel and blotted onto Nytran 
membranes (GE Healthcare). Northern blots were hybridized 
with radioactively labeled probes for open reading frames (ORFs) 
(Prime-IT RmT Random Primer Labeling Kit, Stratagene) or with 
a 5'-phosphorylated SL oligonucleotide CZ4490 CAATATAGTA 
CAGAAACTGTTCTAATAATAGCGTTAGTT. Signals were mea- 



sured using a phosphorimager. Either the signal of the signal recog- 
nition particle (7SL) or the rRNA bands were used for normalization, 
with similar results. 

For high-throughput sequencing, rRNA was depleted using the 
RiboMinus Eukaryote Kit for RNA-seq (Invitrogen). After RNA frag- 
mentation, libraries were constructed and sequenced at the Heidel- 
berg University Bioquant sequencing facility using Illumina kits. 

Bioinformatics 

The analyses were done as described previously (Manful et al. 201 1). 
Seventy-six- nucleotide reads were aligned to the T. brucei TREU 927 
genome, supplemented by all available VSG ORFs from T. brucei 
Lister 427 (Hertz-Fowler et al. 2008) using Bowtie (Langmead 
et al. 2009) and allowing two mismatches. Read counts for each cho- 
sen region were calculated using SAMtools (Li et al. 2009). To avoid 
counting multicopy sequences more than once, we then took a non- 
redundant gene list (Siegel et al. 2010) supplemented by a few addi- 
tional unique genes before calculating the read counts per million 
reads (RPM) and normalizing. Data were normalized from the 
average mRNA SL signal on Northern blots, except in the case of 
the RRP45 RNAi, for which the a tubulin mRNA level was used 
(for explanation, see Results). Half-lives were calculated assuming 
exponential decay, using the formula: 

tl/2 = ln(2) x 30/ - ln("fraction of RPM left at 30 min"). 

In the Tables, the maximum half-life is set to 240 min because of 
the limitations of using single time point. Supplemental Table S2 
(half-lives) contains all annotated genes (coding and noncoding) 
with the exception of rRNAs. In order to analyze the read densities 
across mRNAs, we assigned predominant splice and poly(A) sites 
using all our available data (Manful et al. 2011; Droll et al. 2013). 
Custom PERL and R scripts were written as required. 

Supplemental Table S3 (Differential expression) includes all 
known genes minus rRNAs. Transcripts of genes <250 bp in length 
(the maximum insert size selected during RNA-seq library prepara- 
tions) were not expected to be sequenced, but their precursors might 
have been. Since the precise size cutoff is bound to vary between rep- 
licates, we needed to control for random inclusion of transcripts in 
this size range. Genes of <250 bp with a relative standard deviation 
(standard deviation/mean) between replicates at a steady state of 
>0.3 in either wild-type or knockdown cell lines were therefore ex- 
cluded from the tables. Differential expression was analyzed using 
DESeq (Anders and Huber 2010), with a cutoff P-adjusted-value 
<0.01. 

Class enrichment 

Genes were classified into 72 functional classes, based on annota- 
tions and publications, with 4691 (40%) genes remaining unclassi- 
fied (Supplemental Tables S2, S3). Within a list of genes, class 
enrichment was estimated using a PERL module for the right-hand- 
ed Fisher exact test (Banerjee and Pedersen 2003). The unique cod- 
ing gene list was used for this calculation. Class enrichment heat 
maps were plotted in R using loglO of P-values. 

SUPPLEMENTAL MATERIAL 

Supplemental material is available for this article. 
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